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Abstract 

We give a brief overview of recent developments in Sturm-Liouville the- 
ory concerning operators of transmutation (transformation) and spectral 
parameter power series (SPPS) and propose a new method for numerical 
solution of corresponding spectral problems. 

1 Spectral parameter power series (SPPS) 

Let / G C 2 (a,b) n C 1 [a, b] be a complex valued function and f(x) ^ for any 
x G [a, b]. The interval (a, b) is assumed being finite. Let xq G [a, b]. We 
introduce the infinite system of functions {fk}^L Q defined as follows 

jf(x)XW(x), fcodd, 

fceven, [) 

where 

X^(x) = l, X (n \x) =n ^ X^-^is) (/ 2 (s)) ( ~ ir da, n = l,2,... 
and 

l(°)=l, X^ixj^n T X^^ifis))^ 1 ^ 1 ds, n = l,2,.... 



Together with the system of functions ((T|) we define the functions {V^fcl^o usm § 
the "second half of the recursive integrals, 



ipk(x) 



X^(x) 
77 \ L , k odd, 

XW(x) , 

— , fc even. 

/O) 



1 



As we show below the introduced families of functions are closely related to 
the one-dimensional Schrodinger equation of the form u" — qu = \u where q is a 
complex-valued continuous function. Slightly more general families of functions 
can be studied in relation to Sturm-Liouville equations of the form (py 1 )' + qy = 
Xry. Their definition based on a corresponding recursive integration procedure 
is given in [IT] , [T^] , [5] ■ The system (JJ) is closely related to the notion of the 
L-basis introduced and studied in [7]. 

The following result obtained in [TU] (for additional details and simpler proof 
see [TT] and [H]) establishes the relation of the system of functions {</3fc}^ and 
{ipk}kL to the Sturm-Liouville equation. 

Theorem 1 Let q be a continuous complex valued function of an independent 
real variable x 6 [a, b] and X be an arbitrary complex number. Suppose there 
exists a solution f of the equation 

f" -qf = (2) 

on (a, b) such that f £ C 2 (a, b) D C 1 [a, b] and f(x) ^ for any x G [a, b] . Then 
the general solution u 6 C 2 (a, b) n C^a, b] of the equation 

u" — qu = Xu (3) 

on (a, b) has the form u = c\U\ + CiUi where c\ and C2 are arbitrary constants, 

°° \k 00 yk 

^^M^ and U2 = ^j2kTTy.^ k+1 (4) 

fc=0 v ; fe=0 v ; 

and both series converge uniformly on [a, b] together with the series of the first 
derivatives which have the form 



u'i = f + Y1 (2fe)I [~f lf ' 2k+2k ^ 2k - 1 J and 

00 X k /f \ 

^ = gp^(W 1 + (2 * +1) H' 

The series of the second derivatives converge uniformly on any segment [ai, bj] C 
(a,b). 

The representation (j4|) offers the linearly independent solutions of ([3|) in the 
form of spectral parameter power series (SPPS). The way of how the expansion 
coefficients in (|4]) are calculated via the recursive integration is relatively simple 
and straightforward, this is why the estimation of the rate of convergence of the 
series (@} presents no difficulty, see [12]. Moreover, in [2] a discrete analogue of 
Theorem [1] was established and the discrete analogues of the series (jl]) resulted 
to be finite sums. 
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Remark 2 It is easy to see that the solutions u\ and u 2 defined by satisfy 
the following initial conditions 

ui(x ) = f(x ), u[(x ) = f(x ), 

u 2 (x ) = 0, u' 2 (x ) = l/f(x ). 

Remark 3 It is worth mentioning that in the regular case the existence and 
construction of the required f presents no difficulty. Let q be real valued and 
continuous on [a, b] . Then (0) possesses two linearly independent regular so- 
lutions V\ and v 2 whose zeros alternate. Thus one may choose f = v\ + iv 2 . 
Moreover, for the construction of v\ and v 2 in fact the same SPPS method may 
be used fWj . 

The SPPS representation Q for solutions of the Sturm-Liouville equation (J3J 
is very convenient for writing down the dispersion (characteristic) relations in an 
analytical form. This fact was used in [T] , [5] , [5] , [5] , , [TS] for approximating 
solutions of different spectral and scattering problems. 

2 Transmutation operators 

Let £ be a linear topological space and E\ its linear subspace (not necessarily 
closed). Let A and B be linear operators: E\ — > E. 

Definition 4 A linear invertible operator T defined on the whole E such that 
Ei is invariant under the action of T is called a transmutation operator for the 
pair of operators A and B if it fulfills the following two conditions. 

1. Both the operator T and its inverse T^ 1 are continuous in E; 

2. The following operator equality is valid 

AT = TB (5) 

or which is the same 

A = TBT~ X . 

Very often in literature the transmutation operators are called the transfor- 
mation operators. Here we keep ourselves to the original term coined by Delsarte 
and Lions [5]. Our main interest concerns the situation when A — —-^i + q(x), 

B = —-^pz, and q is a continuous complex- valued function. Hence for our pur- 
poses it will be sufficient to consider the functional space E — C[a, b] with the 
topology of uniform convergence and its subspace Ei consisting of functions from 
C 2 [a, b] . One of the possibilities to introduce a transmutation operator on E 
(see, e.g., [T7]) consists in constructing a Volterra integral operator correspond- 
ing to a midpoint of the segment of interest. Then it is convenient to consider 
a symmetric segment [—a, a] and hence the functional space E = C[— a, a]. It is 
worth mentioning that other well known ways to construct the transmutation 
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operators (see, e.g., [TB], [10]) imply imposing initial conditions on the func- 
tions and consequently lead to transmutation operators satisfying ([5]) only on 
subclasses of E\. 

Thus, we consider the space E — C[— a, a] and an operator of transmutation 
for the defined above A and B can be realized in the form (see, e.g., [TO] and 
[T7]) of a Volterra integral operator 

Tu(x)=u(x)+ I K(x,t)u(t)dt 

J —x 

where the kernel K(x, t) is defined as a solution of a Goursat problem 17, Chap- 
ter 1]. If the potential q is n times continuously differentiable, the kernel K (x, t) 
is n + 1 times continuously differentiable with respect to both independent vari- 
ables. 

In [3] a parametrized family of operators T\, h G C was introduced, given 
by the integral expression 

T h u(x) = u(x) + ( K(x,t;h)u(t)dt (6) 

J —x 



where 



K(x, t;h) = ~ + K(x, t) + ^ [ (K(x, s) - K(x, -s)) ds 
2 2 J t 



(7) 



The operator T; 4 maps a solution v of an equation v" + cu' 2 v = 0, where uj is 
a complex number, into a solution u of the equation u" — q(x)u + lu 2 u = with 
the following correspondence of the initial values 

u(0) = v(0), u'(0) = v'(0) + hv(0). (8) 

Theorem 5 f 1J$ Suppose the potential q satisfies either of the following two 
conditions: 1) q G C 1 [— a, a]; 2) q 6 C[— a, a] and there exists a particular 
complex- valued solution g of non-vanishing on [—a, a]. Then the operator 
Th given by ^ satisfies the equality 
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d 

^2+<?W ) T h [u} = T h 



d 2 

dx 2 ^ 



for any u £ C 2 [— a, a] 



Suppose now that a function / is a solution of @, non- vanishing on [~a, a] 
and normalized as /(0) = 1. Let h := f'(0) be some complex constant. Define 
as before the system of functions {y/c}^ by this function / and by (P). The 
following theorem states that the operator T^ transmutes powers of x into the 
functions ifk- 

Theorem 6 JJJ' Let q be a continuous complex valued function of an indepen- 
dent real variable x € [—a, a], and f be a particular solution of ([2]) such that 
f G C 2 (— a, a) together with l/f are bounded on [—a, a] and normalized as 
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/(O) = 1, and let h :— f'(0), where h is a complex number. Then the operator 
([6]) with the kernel defined by ([7]) (with this particular h) transforms x k into 
ifk(x) for any k <E N . 

Thus, the system of functions {tpk} may be obtained as the result of the 
Volterra integral operator acting on powers of the independent variable. As 
was mentioned before, this offers an algorithm for transmuting functions in 
situations when the explicit form of K(cc, t; h) is unknown. Moreover, properties 
of the Volterra integral operator such as boundedness and bounded invertibility 
in many functional spaces gives us a tool to prove the completeness of the system 
of function {tpk} in various situations. For example, the system {fk)kLo is 
complete in C[— a, a]. 



3 A new method for solving spectral problems 

Theorem[6]together with ((U suggests the following approach for solving spectral 
problems for the Sturm-Liouville equation ([3]) . Here as an example we consider 
the Sturm-Liouville problem with the following boundary conditions though 
obviously the method can be applied in a much more general situation. Thus, 
consider the equation 

u" -qu = -~/3 2 u (9) 

(— (3 2 = A) with the conditions 

u(0) = tt(l) = 0. (10) 

It is well known that an almost optimal uniform approximation of the functions 
sine and cosine on the segment [—1, 1] by polynomials is achieved using the 
Chebyshev polynomials. The corresponding representations have the form [191 
p. 104] 

OO 

smf3x = 2 (-l) m J 2 m+i(/3)IWi(x), (11) 

oo 

COSPX = J {p) + 2j2 {-l) m J2m{P)T2m{x) (12) 
m=0 

where J n denotes the Bessel function of the first kind and T n is the Chebyshev 
polynomial of n-th order, T n (x) = cos(narccosa;). 

From © we have that u(x) — (sin fix) is a solution of © satisfying the 
first of the conditions ([10]). Thus, ff 1 is an eigenvalue of ©, dTUJ) iff T h (sin fix) 
at the point x = 1 equals zero. Notice that 

oo 

T h (sm/3x) =2 (-l) m J2 m +i(P)T h (T 2m+1 (x)) 
-2 J2 (-l) m J2m+i(P)Th {T 2m+l {x)) . 

rn— 
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In the last expression it is not difficult to calculate Th (T 2m +i(x)) with a very 
good accuracy applying theorem[6] Every T\ (T2 m +i(x)) is a linear combination 
of the functions ip^, k = 0, . . . 2m + f and the coefficient corresponding to (pk 
is the same as the coefficient corresponding to x k . Thus, computation of the 
images of the Chebyshev polynomials under the action of the transmutation 
operator Th does not represent any difficulty. Taking into account a very fast 
convergence of the series in (If 2[) one can obtain a very good accuracy in the 
approximation of (sin (3x) with a relatively small N in f) 13|) . 

Thus, the proposed method for solving ©, ([TUf consists of the following 
steps: 1) compute a solution / of @ satisfying the conditions of Theorem[TJ 2) 
compute N + 1 functions ipu] 3) taking the coefficients of x k , k = 0, . . . 2m + 1 
from the expression of T-2m+i compute (T 2m+ i) evaluated at x = 1; 4) find 
zeros of £™ =0 (-l) m ^m+iO^T,, (T 2m+1 ) | x=1 . 

Example g(a^) = e 21 , u(0) = u(ir) — 0. First, making an obvious change 
of variables we reduce the problem to the interval (0, 1) and then compare our 
results with the eigenvalues from [TSJ Appendix A]. Taking N = 18 we obtain 
the results reported in the following table. The recursive integration involved in 
the construction of the functions tpk required in (| 13|) was performed by means of 
the Spline Toolbox of Matlab. On each step the integrand was converted into a 
spline and then integrated using the command fnint. A surprising accuracy with 
a very small number of functions participating in (If 31) was achieved in all the 
considered numerical tests. Here due to the restrictions in space we illustrate this 
with the only example in which N = 18. A phenomenon which we observed but 
up to now have not found an explanation is a surprisingly good approximation 
of higher order eigenvalues as can be seen from the table. Apparently the 
expression (fT3| not only approximates well the solution of the equation for 
relatively small values of /3 but also evaluated at x = f contains information on 
the asymptotics of the eigenvalues of the problem. It is worth mentioning that 
on a usual netbook computer the computation of the eigenvalues (see the table 
below) implemented in the Matlab takes only several seconds. 
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Order of the eigenvalue 


Eigenvalues from 1181 


Computed eigenvalues 


1 


4.8966693800 


4.896669377 


2 


10.045189893 


10.045189883 


3 


16.019267250 


16.019267268 


4 


23.266270940 


23.266270969 


5 


32.26370704 


32.26370710 


6 


43.2200196 


43.2200184 


7 


56.18159 


56.18161 


8 


71.15299 


71.15255 


9 


88.1321 


88.1386 


10 


107.11 


107.05 


11 


128.10 


128.48 


17 


296.07 


296.52 


28 


791.05 


790.99 


43 


1856.05 


1856.53 


50 


2507 


2500 
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